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ABSTRACT 

A  numerical/analytic.al  procedure  has  been 
developed  which  yields  the  optimum  amount  of  rein¬ 
forcement  for  a  given  hole  shape  in  a  large  elastic 
plate  under  prescribed  boundary  stresses.  This  pro¬ 
cedure  is  based  on  determining  the  usual  two  complex 
potentials  which  describe  the  entire  stress  field, 
constructing  the  strain  energy  density  function  in 
terms  of  the  unknown  amount  of  reinforcement, 
integrating  this  function  around  the  opening  boundary, 
and  finally  minimizing  this  integral  with  respect  to 
the  reinforcement. 

The  method  is  first  developed  for  a  general  hole 
shape  and  then  demonstrated  in  some  detail  for  a  cir¬ 
cular  and  a  square-like  opening. 


INTRODUCTION 


The  use  of  reinforcements  at  the  boundaries  of  openings  is  standard  prac¬ 
tice  in  ship  and.  aircraft  construction.  The  amount  of  reinforcement  to  be 
used  is  determined,  in  most  cases,  by  somewhat  arbitrary  decisions  regarding 
the  percentage  reduction  it  produces  in  the  boundary  stress  maximums. 
Although  intrinsically  there  is  nothing  wrong  with  this  approach,  it  seems  to 
lack  a  rational  basis  or  criterion  for  helping  the  structural  designer  select 
a  specific  amount  of  reinforcement.  The  procedure  presented  in  this  paper  is 
an  effort  to  provide  such  a  basis.  Rather  than  determining  the  stress  field 
corresponding  to  a  given  geometry  of  the  opening  and  of  the  reinforcement, 
this  procedure  seeks  to  determine  that  reinforcement  which  minimizes  a  certain 

meaningful  integral  related  to  the  boundary  strain  energy.  In  this  way  it  is 

1-3 

an  inverse  elasticity  problem.  Some  investigators  ,  in  studying  the 
related  problem  of  optimizing  unreinforced  notch  shapes  in  plates,  have  con¬ 
cluded  that  uniform  tangential  stress  at  the  notch  boundary  would,  in  general, 
lead  to  the  smallest  stress  concentrations.  Intuitively,  it  appears  justified 
to  assume  that  in  the  case  of  reinforced  notch  boundaries  the  requirement  of 
uniformity  of  boundary  stresses  and/or  stress  related  quantities,  such  as 
strain  energy  etc.,  at  the  notch  boundary  would  lead  to  more  desirable 
designs.  The  optimization  rationale  used  in  the  present  procedure  is  based  on 
this  argument . 


I 

It 
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In  general,  stress  analysis  of  a  non-circular  opening  reinforced  with  a 
thin  member  of  uniform  cross  section  is  very  difficult  because  it  requires  the 
satisfaction  of  a  boundary  condition  which  contains  an  irrational  term.  How¬ 
ever,  the  use  of  MACSYMA  (a  symbolic  manipulation  language  developed  at  MIT 
and  in  use  at  DTNSRDC)  makes  it  possible  to  solve  such  problems,  since  a 
larger  number  of  terms  can  be  retained  and  manipulated  in  various  expansions 
without  losing  track  of  them  in  the  enormously  long  and  complex  algebraic 
expressions. 

The  general  form  of  the  boundary  condition  used  in  this  method  is 
developed  in  reference  4.  The  special  cases  of  reinforcement  of  a  circular 
opening  and  a  square-like  opening  are  discussed  in  greater  detail  and  actual 
numerical  results  are  included. 

This  paper  is  a  logical  extension  of  the  work  described  in  reference 
which  dealt  with  the  optimization  of  the  shape  of  a  class  of  unreiniorced 
openings  in  large  plates. 


MATHEMATICAL  PRELIMINARIES 

An  opening  of  general  shape  in  a  large  elastic  plate  of  isotropic 
material  of  unit  thickness  is  reinforced  by  a  thin  member  of  cross-sectional 
area  A  capable  of  withstanding  axial  forces  only.  Then,  if  the  opening  can  be 
mapped  into  a  unit  circle  by  a  function  z(0,  the  equivalence  of  complex 
forces  between  those  in  the  reinforcement  and  those  in  the  plate  at  the  open¬ 
ing  boundary  (4  =  0)  is  given  by  4 

i .(a)  +  4>'(a)  +  =  a  P  \l  Z-'--q)  +  c  O) 

z ' (0)  '  z ' (o) 

where  4>  (  o  )  and  P  (  o  )  are  the  values  of  the  functions  +  (  C  )  and 

ip(r)  at  the  opening  boundary,  P  is  the  axial  force  in  the  reinforcement, 

,  iB 

o  equals  e  ,  ana  c  is  an  arbitrary  constant.  Equivalence  of  tangential 
strain  at  the  boundary  requires 

P  -  A  (o  -  ye  >  (2) 

p  a 
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t 


and  the  elastic  equilibrium  of  a  boundary  element  requires 


a  =  P 
a 


I  z  z 


(3) 
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1 

J“7=i 
\  z  z 


SP 
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where  the  argument  o  of  the  function  z  has  been  omitted  for  brevity  and 
prime  indicates  differentiation  with  respect  to  a 

The  well-known  relations  between  the  stresses,  o  ,  o„  ,  and  t  .  and 

*  a  *  B  aB 

the  functions,  4>  (  C  )  and  ij  <  (  t),  and  Equations  (2)  and  (3)  can  be  used  to 

show  that 


Vz'  z  '  j^A  (1+v)  +  Vz'  z  '  J  =  2A(ij),z'  +  4> '  z  * ) 


(4) 


Equation  (1)  can  be  modified  to 


(4>  +  iji)  z’  +  z  $ 1  =  o  P  Vz 'z 


(5) 


where  the  argument  a  has  been  omitted,  for  brevity,  for  the  $  ,  ib  ,  and 
z  functions;  prime  indicates  differentiation  with  respect  to  the  argument  a  ; 
and  a  bar  (-)  represents  the  complex  conjugate  of  the  function.  Since  P  is  a 
real  quantity,  it  is  possible  to  represent 


Gr  ■  I 

n=0 


e  (o 
n 


(6) 


where  the  c  are  real.  Functions 
n 


following  form: 


>  (  C )  and  ♦  (  C )  are  known  to  have  the 

♦  (c)  =  Sc  +  X 


(7) 


♦  (c)  -  DC  +  Y-r 
cn 


y  t 


•*  ' 


where 


p  +  q 

S  =  - 

4 
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and  p  and  q  are  uniform  stresses  at  infinity  at  an  angle  9  to  the  x-  and  y- 
axes,  respectively.  The  mapping  function  z(  4 )  can  be  conveniently 
represented  by 


z(?) 


(8) 


wnere  m  are  known  constants  and  R  is  a  size  factor  which  will  be  assumed 
n 

unity. 


The  problem  can  now  be  stated  in  two  parts:  (a)  For  a  given  mapping  func¬ 
tion  z(  4  ) ,  determine  the  functions  <J>  (  4  )  and  ’l'  (  4  )  as  a  function  of  the 
area  of  reinforcement  A;  and  (b)  determine  that  value  of  A  which  minimizes  a 
certain  meaningful  integral  related  to  the  boundary  strain  energy. 

Clearly  the  solution  of  part  (a),  which  forms  the  basis  for  the  solution 
of  part  (b),  would  require  the  determination  of  an  and  b^  such  that  Equa¬ 
tions  (4)  and  (5)  are  satisfied. 

Substitution  of  Equations  (6),  (7)  and  (8)  into  the  boundary  condition. 
Equation  (5),  and  its  ensuing  solution  generates  one  set  of  linear  simultane¬ 
ous  Equations  between  a„  and  c  and  another  set  between  a  b 

n  n  n  n 

and  cn  ,  both  of  which  satisfy  this  boundary  condition.  Values  of  a^  in 

terms  of  cn  can  be  determined  from  the  first  set  and  substituted  in  Equation 

(4>  to  determine  cn  in  terms  of  A  by  using  an  expansion  of  the  irrational 

term  Vz  'z '  .  Taylor  series  expansion  of  this  term  at  C  =  u  appeared  to 

be  a  reasonable  choice,  and  it  was  seen  that,  for  the  demonstration  problem  of 

6 

a  square-like  hole ,  an  accuracy  of  1  in  10  could  be  obtained  at  the  hole 
boundary  by  using  the  first  eight  terms  after  the  constant  term.  Once 
the  cn  are  known,  the  two  sets  of  simultaneous  equations  can,  in  principle, 
be  solved  to  obtain  an  and  bn  in  terms  of  the  unknown  quantity  A.  The 


4 


basic  information  for  the  optimization  problem  is  then  available 


OPTIMIZATION 


The  strain  energy  density,  Vq  ,  of  a  plane  elastic  system  is  given  by 


II 

o 

> 

1 

2E 

r  2  2 

a  +  a.  -  2va  a„ 

6  a  B 

+  20+V>Tag2] 

(8) 

With  Equations 

(2) 

and  (3),  this  expression 

for  V  can  be  transformed  to 

0 

II 

0 

> 

1 

4E 

te +  p2 

A  z  z  J 

,  4(l+v)  /  8P\2 

z'z’  W/ 

(9) 

The  area  of  reinforcement,  A,  as  it  occurs  in  Equation  (9;  is  a  dimen¬ 
sionless  quantity  since  it  has  been  divided  by  the  unit  thickness  of  the  plate 
as  well  as  by  the  unit  size  factor  K  of  the  mapping  function. 

For  specific  loading  conditions,  the  strain  energy  density,  V  ,  a  func- 

o 

tion  of  B  and  A  only,  can  be  integrated  with  respect  to  8  from  U  to 
2  it  to  obtain  an  integral,  I,  which  can  be  interpreted  as  the  strain  energy 
xn  a  thin  region  of  variable  thickness  around  the  opening.  This  integral  can 
now  be  minimized  with  respect  to  A  to  obtain  that  value  of  A  which  leads  to 
the  "most  uniform"  distribution  of  Vq  around  the  opening.  The  term  "most 
uniform"  conveys  the  meaning  that,  of  all  the  possible  parametric  variations 

of  V  corresponding  to  a,  this  particular  distribution  would  have  strongly 

°  5 

attenuated  peaks. 


A  distinction  must  be  made  between  this  integral,  i,  and  the  actual 
strain  energy,  V,  in  a  thin  slice  around  the  opening.  The  strain  energy 
itself  will  be  obtained  by  integrating  Vq  with  respect  to  ds.  Since  ds  is 
given  by  V  z '  z '  dB  ,  the  strain  energy  will  be 


(10) 


i 
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It  can  be  expected  that  a  minimization  of  V  with  respect  to  A  will  lead 

to  the  "most  uniform"  distribution  of  V  V z ' z  ’  rather  than  that  of 

o 

around  the  opening.  The  integral,  1,  was  previously  calculated,  within  a 
constant,  in  reference  (5)  for  a  square-like  unreinforced  opening  to  determine 
its  optimum  shape.  For  S  =  .25  and  D  =  -.5  the  integral , I  ,was  found  to  have  a 
minimum  value  of  0.3627  at  m  =  -.05.  This  value  of  m  produced  a  very  desir¬ 
able  stress  distribution  and  the  lowest  stress  concentrations  since  for  an 

2 

unreinf orced  opening  the  integrand,  ,  of  I  was  °g  .  For  comparison, 

when  this  case  was  later  repeated  by  minimizing  V  such  that  the  integrand  was 
■I,-  V z'z'  ,  the  minimum  value  of  V  within  the  same  constant  was  found  to  be 
0.3705  at  m  =  -.07.  The  corresponding  stress  distribution  around  the  opening 
is  less  desirable  than  when  m  =  -.05.  in  fact, the  highest  stress  concentra¬ 
tion  for  m  =  -.07  was  2.505  versus  2.472  for  m  =  -.05.  Never theiess , 

the  o  '’  V z ' z  1  distribution  can  be  expected  to  be  "most  uniform"  at  m  =  -.07 
and  the  J  -  distribution  is  the  "most  uniform"  at  m  =  -.05.  f  rom  an 

O 

engineering  stand  point  the  "most  uniform"  -  distribution  will 

obviously  be  preferable.  Other  more  mathematical  reasons  for  choosing  to 
minimize  the  integral  of  o  are  included  in  reference  5. 

The  preceding  discussion  justifies  the  use  of  the  integral,  1,  rather 

than  the  strain  energy,  V,  as  the  quantity  to  be  minimized  for  determining  the 

optimum  A.  in  summary,  minimization  of  1  should,  in  general,  be  expected  to 

lead  to  the  "most  uniform"  V  -  distribution  and  minimization  of  V  to  the 
r— rr  ° 

"most  uniform”  V  Vz'z'  -  distribution.  It  can  also  be  expected,  as  evidenced 
in  tiie  case  of  the  unreinforceo  square  opening,  that  "most  uniform"  - 

distribution  would  result  in  smaller  V  values. 

o 

At  this  point  it  should  be  mentioned  that  consideration  was  also  given  to 
determining  optimum  A  based  on  attenuation  of  other  relevant  quantities  such 
as  the  distortion  energy  density  tequivalent  to  attenuation  of  the  octahedral 
stress),  the  dilatation  energy  density,  the  maximum  shear  stress,  and  the  max¬ 
imum  principal  stress.  All  these  quantities  except  the  last  are  either  as 
difficult  or,  in  fact,  simpler  to  attenuate  than  the  strain  energy  density. 
Since  the  attenuation  of  the  maximum  principal  stress  is  in  general,  consider¬ 
ably  more  difficult  to  perform  (and  not  necessarily  most  desirable),  t lie 
strain  energy  density  was  selected,  in  the  case  of  the  square-like  hole,  as  a 


6 


typical  quantity  to  be  attenuated.  in  the  case  oi  the  circular  hole  other 
quantities  were  also  investigated  for  comparison. 


EXAMPLES 


L.  Circular  Opening 

In  this  case  the  irrational  terms  in  Equations  t  3 )  ,  (4),  and  (5)  disap¬ 
pear,  so  the  manipulations  are  rather  straight  forward.  Another  point  to  be 
noted  is  that  tot  a  circular  opening  the  strain  energy,  V,  is  the  same  as  Lite 
integral  i.  Omitting  ail  the  intermediate  manipulations,  the  strain  enemy. 
V,  can  be  given  by 

2(1  +  0. 91A2)  S2  +  U  +.1.1^3 1A.V.  (Ha) 

(1  +  1.3A)2  (1  +  3.3A)2 


tie  energy  oi  distortion,  Vn  ,  or  the  integral  of  the  square 

2 

tress,  TQct  ,  can  be  given  by 


ol  *.K  tailed  1  . . 


■I 


2 it  3 . 9t 


2E 


oct  ,a  13.867tr 
do  — 


2E 


2(1-0. 4A+0.79A2)S2  (1-0.4A+12.79A2)D2 

9  +  9 


(l  +  1.3AV 


n  +  1  iai 


(lib) 


t lie  energy  due  to  change  in  volume,  Vv  ,  can  be  given  by 


V 

v 


2. 1 3  3tt 
2E 


(1  +  1 . 3A)  2 
(1  +  3. 3A) 2 


(11c) 


and  finally  the  integral,  T,  of  the  square  of  maximum  snear  stress  can  lx 
given  by 


T  =  16ir 


2(1  -  1.4A  +  0.49A2)S2  + 


(1 


1.4A  +  16.49A2)!)2 


(1  +  1 . 3A) ‘ 


(1  +  3. 3A) ‘ 


(lid) 
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TABLE  1  -  OPTIMUM  VALUES  OF  REINFORCEMENT  FOR  A  CIRCULAR  HOLE 
SUBJECT  TO  VARIOUS  TYPES  OF  LOADS: 


ATTENUATED 

QUANTITY 

STRAIN 

ENERGY 

DENSITY 

DISTORTION 

ENERGY 

DENSITY 

VOLUME 

ENERGY 

DENSITY 

MAX  SHI 
STRES 
SQUARJ 

£AR 

s 

ED 

Load 

C/\StS 

V 

o 

rm 

VD 

A 

V 

V 

■ 

KBi 

A 

isot  ropic 
s=u .  b ;  D=U 

2.8l1 

‘  2E 

1.4286 

1.  733-rr 
2E 

1.4286 

1.067 n 
2E 

B 

o 

1.4286 

i 

Uniaxial 

S=U.25;  D=-U. 5 

3.074-rr 

2E 

.4429 

2.584ir 

2E 

.4033 

0.349ir 

2E 

** 

00 

2. 717H 

.  3422 

Pure  Shear 

S=U ;  U=-i 

8 . 15ir 
2E 

.2918 

7.072ir 

2E 

.2602 

0.331it 

2E 

** 

00 

8tt 

.2127 

*  -  Indicates  that  Vv  is  independent  of  A 

**  -  Infinity  refers  to  a  rigid  reinforcement 


The  quantities  V,  VR  ,  V  >  anc*  ^  were  minimized  to  determine  the 
optimum  values  of  A.  These  computations  were  performed  for  a  number  of  load 
cases  as  shown  in  Tabie  1.  The  values  of  A  at  which  the  minimums  occurred  and 
the  corresponding  minimized  quantities  are  included  in  this  table. 

'Hie  tact  that  such  optimum  reinforcements  can  exist  is  intriguing.  A 
popular  definition  of  A  is  the  ratio  of  area  replaced  as  reinforcement  to  the 
area  removed  by  the  opening.  Thus  in  the  uniaxial  case,  if  approximately 
44. 3%  of  the  area  removed  by  the  opening  is  replaced  by  the  reinforcement,  the 
resulting  stress  distribution  would  be  very  desirable,  i.e.,  one  that  would 
result  in  a  relatively  smooth  boundary  strain  energy  density  distribution. 
The  actual  stress  distributions  corresponding  to  the  cases  in  Table  1  can  be 
obtained  using  equations  derived  in  reference  6. 
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II.  Square  Opening 

The  algebraic  manipulations  in  this  case  are  considerably  more  involved 
and  lengthier  than  for  the  circular  opening.  A  square  opening  described  by 

z  -  a  -  ^5  (12) 


was  chosen  for  this  example  because  of  its  optimum  square-like  shape  .  Solu¬ 
tion  of  Equations  (4)  and  (5)  and  integration  of  Equation  (9)  required  the  use 


of  the  irrational  terms  (z'z') 


(z'z ')  2  ,  and  (z' z' )  ^ 


which 


convenient] 


represented 
4n  ,  -4n, 


Taylor 


type  ,4  G4n  (o  +  a  ),  c.f.  Equation  (15). In  all  cases  truncation 
after  eight  terms  could  maintain  an  accuracy  of  1  in  106  in  the  boundary 
values  of  these  expansions.  Equation  (5)  produced  two  sets  of  equations,  the 
solution  of  one  of  which  is  given  by 


a9  — .<5  .ma  =  c„  +  3ma„  .  +  6  „mS  -  6  D 
2n-l  n,l  1  2n  2n+3  n,2  n,l 


a-  =0 
2n 


;  n  >  1 


where  m  =  -0.05.  it  is  easy  to  see  from  Equation  (13)  that,  if  either 

a  or  c  is  truncated,  the  values  of  individual  a  can  be  easily  deter- 
n  n  n 

mined  in  terms  of  c  .  Because  of  the  small  value  of  m,  the  final  results 


converged  within  1  in  10  when  the  a  beyond  n  =  13  were  truncated,  unce 

n 

the  a  were  determined,  the  t,  could  be  determined  using  the  other  set  ot 
n  n 

equations.  However,  this  procedure  was  not  necessary  since  the  alternative 
procedure  for  accomplishing  the  solution,  described  in  the  next  paragraph,  was 
more  direct. 


The  a  and  Taylor  expansion  of  V z  z '  were  substituted  in  Equation  (4) 
n 

to  obtain  c  ,  given  by 

n  \ 

C2n-1  =  0  ;  n  >  1 


qk  A  ;  4  >  n  >  0 


(i  *k)s/  i’k*k 

\  k*l  /  k*0 

"IS  P(4n_2)'k  Ak)D/  ^  r 

\  k-1  }  k-0 


k  ;  4  >  n  >  1 


where  p  ,  q  ,  and  r  are  known  coefficients  but  are  not  included  here, 
n  n  n 


Once  the  c^  were  obtained,  the  problem  was  basically  solved  as  a  func¬ 
tion  of  the  area  of  reinforcement  A,  The  were  substituted  in  Equation  (6) 

to  obtain  P  vf7  z  which,  in  turn,  was  substituted  in  Equation  (9)  to 
obtain  Vq  .  The  stresses  °a  »  °g  >  anc*  Tag  could  be  obtained  using 

Equations  (2)  and  (3).  In  order  to  integrate  V  ,  the  following  Taylor 

o 

expansions  were  used: 


8 


(z 

,  /  4n 

L4n(° 

-4n. 
+  o  ) 

n=0 

8 

(z 

o'1  -  z 

,  4n 
M4n  ° 

,  -4n. 

+  0  ) 

R=0 

These  expansions,  witli  the  now  known  p  \4  '  z ' ,  yielded 


8 

P  =  Z  P2ncos  2nB 
n=0 

,  -  .  _! 


8 


(z  z  )  =  /  Q2ncos  2nB 


n=0 


(15) 


(16) 


which  were  used  to  obtain  the  integral 

"2  7T  21 


-f 


vn  dB  =  +  2(1-v2)T  +  4(l+v)I 

'0  0  A  1  3 


(17) 


where 


L1  2E  (  Po  +  2 


1  y 

2 


\  ■  w  I  %2  * 


n=l 

8 


2n 


n=l 


2n 


I  =  —  M 
3  2E 


j  i  o»v2  ♦  z 

n=l  n=l 


8 


P  2  + 

4n  2n  '‘4n‘'4n, 

n=l 


z 


(18) 


The  R  ,  given  as  sums  of  products  p  p  ,  are  coefficients  of  cos  n8  in 

n  /sp\2  m  n 

the  expression 


/spy 

V38/ 


.  The  minimum  value  of  1,  which  in  Equation  (17)  is 
a  function  of  A,  S,  and  D  only,  was  determined  by  simply  evaluating  I  at 
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various  values  of  A  for  a  selected  set  of  values  of  S  and  D.  In  fact,  the 


three  cases  considered  were  the  isotropic  case  (S  =  .5,  D  =  0) ,  the  uniaxial 

case  (S  *  .25,  D  =  -.5),  and  the  pure  shear  case  (S  =  0,  D  =  -1).  The 

corresponding  A's  at  which  the  minimum  values  of  I  occurred  were  0.90,  0.41, 

and  0.29.  These  minimum  values  of  I  could  be  directly  compared  to  those  of 

2E 

reference  (5),  if  they  were  first  multiplied  by  .  The  I  values  for 

these  three  cases  were  found  to  be  .2453,  .1920,  and  .4771,  respectively.  It 
should  be  noted  that  since  the  value  of  m  for  all  three  cases  was  -0.05,  the 
initially  unreinforced  shape  is  optimized  only  for  the  uniaxial  case ^  .  The 
minimization  of  I  with  respect  to  A  for  ail  three  cases  demonstrates  that 
there  is  also  an  optimum  amount  of  reinforcement  for  those  openings  which  do 
not  initially  have  an  optimum  shape.  Table  2  gives  the  minimized  values  of  1, 
within  the  constant,  and  the  corresponding  values  of  A  for  the  uniaxial  case 
only . 


TABLE  2  -  MINIMIZED  VALUES  OF  I  AT  OPTIMUM 
VALUES  OF  A  FOR  A  SQUARE  OPENING 


A 

4.JC, 

16tt  ’ 

0.000001 

.36268 

0.1 

.24854 

0.2 

.20936 

0.3 

.19546 

0.4 

.19199 

0.412 

.19195 

0.5 

.19324 

0.6 

.19671 

0.7 

.20120 

0.8 

.20610 

0.9 

.21109 

1.0 

.21600 

the  value  of 

2E 

16tt 

I  correspond!! 

same  as  that  of  reference  (5)  for  A  »  0  (unreinforced). 


.00001  *  0  is  the 
Table  2  also  snows 


the  variation  of  I  with  respect  to  A  and  its  minimum  point  at  A  =  .412. 
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Figure  1  shows  the  stresses  a 
like  opening  (m  =  -.05)  with  A 


and  t  „  for  a  reinforced  square- 
a8 


a  '  8 

.412  and  also,  for  comparison,  the  equivalent 

and  the  value 


stress,  o  ,  based  on  the  strain  energy  theory  of  failure 
E 

of  a  for  the  unreinforced  opening  subjected  to  uniaxial  loading.  The  maximum 

8 

value,  1.43,  of  a  in  the  reinforced  case  occurs  approximately  at  6  =  80°and 
p 

in  the  unreinforced  case  this  maximum  2.47,  occurs  at  approximately  8  =  70°. 


Since  the  amount  of  area  removed  is  not  a  constant  by  definition,  in  the 
case  of  a  square  opening  with  rounded  corners,  the  percentage  of  area  replaced 
can  not  be  a  constant  either.  For  convenience,  however,  it  can  be  referred  to 
the  half-width  of  the  square.  Thus  the  value  of  A  should  be  divided  by  1  +  m, 

i.e.,  0.95,  and  multiplied  by  100  to  obtain  the  reinforcement  expressed  as  a 
percentage  of  the  area  removed. 


CONCLUSIONS 

1.  It  has  been  shown  that  for  openings  in  plates  an  optimum  amount  of  rein¬ 
forcement  exists  which  corresponds  to  a  minimum  value  of  a  strain  energy 
related  integral. 

2.  It  can  be  conjectured  that  the  values  of  strain  energy  densities  around  a 
reinforced  opening  are  relatively  smooth  when  this  strain  energy  related 
integral  is  a  minimum. 

3.  The  actual  optimum  values  of  reinforcement  for  a  circular  opening  were 
found  to  be  142%,  44.3%,  and  29.2%  of  the  area  removed  by  the  opening  for 
isotropic,  uniaxial,  and  pure  shear  loadings,  respectively. 

4.  The  optimum  values  of  reinforcement  for  a  square-like  opening  were  found 
to  be  94.7%,  43.4%,  and  30.5%  of  the  area  removed  by  the  opening  for  iso¬ 
tropic,  uniaxial,  and  pure  shear  loadings,  respectively. 
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FIGURE  1  -  Various  Boundary  Stresses  Around  a  S<juare-Like  Opening 

With  Optimum  Reinforcement 
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